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Abstract 

A dynamical approach is developed to predict the vrA^ scattering amplitudes 
starting with the constituent quark models. The first step is to apply a vari- 
ational method to solve the three-quark bound state problem. The resulting 
wave functions are used to calculate the A^* ttN, r]N, ttA vertex functions 
by assuming that the vr and r] mesons couple directly to quarks. These vertex 
functions and the predicted baryon bare masses then define a Hamiltonian for 
ttN reactions. We apply a unitary transformation method to derive from the 
constructed Hamiltonian a multi-channel and multi-resonance reaction model 
for predicting the irN scattering amplitudes up to = 2 GeV. With the pa- 
rameters constrained by the A(1232) excitation, we have examined the extent 
to which the vrA^ scattering in 5ii channel can be described by the constituent 
quark models based on the one-gluon-exchange or one-meson-exchange mech- 
anisms. It is found that the data seem to favor the spin-spin interaction due 
to one-meson-exchange and the tensor interaction due to one-gluon-exchange. 
A phenomenological quark-quark potential has been constructed to reproduce 
the 5ii amplitude. 



PACS Numbers: ll.80.Gw, 12.39.Jh, 12.40.Yx, 13.75.Gx, 14.20.Gk 
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I. INTRODUCTION 



The constituent quark models have long been used to investigate the structure of nucleon 
resonances. Most of the earlier works were based on phenomenological forms of residual 
quark-quark {qq) interactions. With the development of Quantum Chromodynamics(QCD), 
a more fundamental approach was developed by assuming that the residual qq inter- 
actions can be parameterized as the Fermi-Breit form of one-gluon-exchange mechanism . 
In recent years, an alternative approach has been developed p|-p^ based on the assumption 
that the residual gg-interaction is due to the exchange of octet Goldstone bosons. With 
appropriate phenomenological tuning, both approaches can reproduce the general structure 



of the baryon spectra listed by Particle Data Group (PDG) [0. Some attempts [|T^-[T5 



have also been made to develop hybrid models including both one-gluon-exchange and one- 
meson-exchange quark-quark interactions. To make progress, it is important to develop an 
approach to distinguish all of these constituent quark models using additional experimental 
data. 

The main question we want to address in this work is how the constituent quark models 
can be tested against the vriV scattering data. It is common to compare their predicted 
masses and decay widths with the data listed by PDG. All calculations of decay widths have 
been done |T^-^ perturbatively. For example, the width of the decay of N* into a vriV state 
is calculated by evaluating the matrix element (A^*|0(A;)|A^) with an appropriate operator 
0{k) describing how pions are coupled to quarks. The interactions between the outgoing 
mesons and baryons are neglected. It has been found that such a perturbative calculation 
can at best describe the general qualitative trend of the data, but not the quantitative 
details. For example, the widths of A^(1535) — > 7[N,7]N,7rA predicted by Ref. |20| are 
(14.7 ± 0.5, 14.6 ± 0.4, 1.4 ± 0.3) MeV^/^ ^ ^i^j^i^ 

seem in quantitative agreement with 
the empirical values (8.0 ± 2.8, 8.1 ± 0.8, 0.) MeV^/^ determined in Ref. [ggf. 

It is important to note here that the PDG's values are extracted from the experimental 
ttN amplitudes which contain both resonant and non-resonant components. In most partial 
waves, the non-resonant mechanisms are important ; one can see this from the fact that most 
of the resonances identified by PDG are in fact not visible in vriV and cross section data. 
By the unitarity condition, therefore the extracted resonance parameters "inherently" con- 
tain non-resonant contributions. Furthermore, the separation of non-resonant components 
from the full amplitudes is a model-dependent procedure. The available amplitude analyses 
2|-|3 have yielded very different resonance parameters in many cases. Clearly, except in 



a region where the non- resonant contributions are negligibly small, the comparison of the 
PDG values (or values from other amplitude analyses) with the decay widths calculated 
perturbatively from the constituent quark models could be very misleading. In particular, a 
perturbative calculation of decay widths is obviously not valid for cases in which two nearby 
resonances in the same partial wave can couple with each other through their coupling with 
the meson-nucleon continuum. Similar precautions must also be taken in comparing the 
predicted masses with the PDG values. These issues concerning the comparison of the PDG 
data with the predictions from constituent quark models were discussed in Ref. [p8| . 

To have a more direct test of constituent quark models, we will explore in this work a 
nonperturbative approach that takes account of the unitarity condition and can relate the 
TiN scattering amplitudes directly to the predicted internal quark wave functions of baryons. 
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Our approach is guided by a dynamical model of nN and reactions developed in Ref. 
29| (SL model). It was shown there that the vrA^ and reactions up to the A(1232) 



energy region can be described by the following Hamiltonian 

H = Hq + TA^TTN,'yN + VttN,ttN + VTTN,'yN , (l) 

where rA^7rAr,7Af describes the A ^ ttN, •jN transitions, and Va,i3 are the non-resonant in- 
teractions. It was found that the unitarity condition and the non-resonant interaction Va,f3 
can shift the mass of A by about 60 MeV and account for as much as about 40 % of the 
Ml strength of the A 'jN decay. This provides an explanation of a long-standing dis- 
crepancy between the Ml value predicted by the constituent quark model and the PDG 
value. It is therefore natural to conjecture that Ho of Eq. (Q) can be identified with the 
model Hamiltonian of a constituent quark model, and rA-*7rAf,7Ar correspond to the decay 
amplitudes calculated perturbatively using the resulting baryon wave functions. This as- 
sumption is then similar to what was used in a dynamical study [BO] of the vrA^ scattering 



amplitude in channel within a constituent quark model. However, Ref. did not 
consider the non-resonant interaction Va,f3 and employed very simple internal wave functions 
for the A^*(S'ii) states. In this work, we will extend the Hamiltonian Eq. (|l]) to consider 
the multi-channel and multi-resonance cases. As a first step, we will focus on the Su chan- 
nel. We will concentrate on analyzing the dynamical content of our approach in this rather 
complex channel. 

There have been some attempts to understand the constituent quark models within QCD. 
Manohar and Georgi ||^ argued that in the kinematic region between the chiral symmetry 



breaking scale A^sb ~ 1 GeV and the QCD confinement scale Aqcd ~ 0.1 — 0.3 GeV, the 
effective theory for hadrons is defined by the Lagrangian 

L{g, G, 0) = ^Pizd^-j, - m,)^ - igt/jG^j^ij - ||V'7^5V'5m</> + ■ ■ ■ ■ (2) 

where is the constituent quark mass, G^i and are the fields for constituent quarks , 
gluons, and Goldstone bosons, respectively. The most crucial dynamical assumption of this 
approach is that the Goldstone bosons are coupled directly to the constituent quarks by the 
flavor SU(3) symmetry characterized by the coupling constant 5f^/(2/). This is consistent 
with the notion that the Goldstone bosons result from the spontaneously breaking of the 
approximate chiral symmetry that characterizes the QCD Lagrangian. 

The Lagrangian Eq. implies that the constituent quarks could interact with each 
other through both the exchanges of gluons and Goldstone bosons. The usual one-gluon- 
exchange (OGE) and one-meson-exchange (OME) gg-interaction can be calculated from 
Eq. (I) using perturbation theory. This conjecture is supported by a recent Lattice QCD 
calculation [Q. In Ref. |32|, it was found that the mass splitting between N and A is 



largely due to the meson-exchange mechanism. The need of both flavor-independent (i.e. 
OGE) and flavor-dependent (i.e. OME) gg-interactions is also suggested by the baryon mass 
formula determined in an algebraic approach In this work, we will take this point of 



view and will consider a general constituent quark model which can have both OGE and 
OME mechanisms. This model can be reduced to the previously developed OGE or OME 
models in some limits. 

The meson-quark couplings were also included in other hadron models such as the chi- 
ral/cloudy bag models |3^. However the Lagrangian Eq. (^ with ~ 200 — 300 MeV for 



3 



up and down quarks is closest to a framework within which one can hope to understand the 
dynamical origins of the most often used constituent quark models based on gg-potentials. 
Perhaps one can apply the dynamical approach developed in this work to also examine other 
hadron models. But this is beyond the scope of this paper. 

In section II, we present the model Hamiltonian for baryons and our method for solv- 
ing the three-quark bound state problem. The calculations of meson-baryon-baryon vertex 
functions are given in section III. A formalism for vrA^ reactions is developed in section IV. 
In section V, we analyze the dynamical content of our approach within a simple model. The 
results and discussions are presented in section VI. Section VII is devoted to conclusions 
and discussions on future developments. 



II. INTERNAL STRUCTURE OF BARYONS 

A. Baryon Hamiltonian 

In this exploratory study, we take the simplest, but often used, approach and assume that 
the baryon structure can be described in terms of three nonrelativistic constituent quarks. 
The model Hamiltonian is of the following familiar form 

hB=K + Konf + Vqq . (3) 

The kinetic energy K is defined by 

where Pi is the momentum of the ith quark and P is the center of mass momentum of 
the three-quark system. The constituent quark mass rriq is taken to be 340 MeV, which is 
close to the value used for describing the nucleon magnetic moments within the simple (Os)^ 
configuration. The confinement potential Konf is assumed to be of the usual linear form 

Konf = "c'^iJ ' (^) 
i<j 

where Vij = \ri — rj\ and ttc is a constant. 

For the residual gg- interaction Vqq in Eq. (|^), we first consider the cases that either the 
one-gluon-exchange(OGE) model or the one-meson-exchange(OME) model is used. Both 
models are derived from taking the static limits of the one-particle-exchange Feynman am- 
plitudes. For the OGE model, following the previous works, we drop its spin-orbit component 
and retain only spin-spin and tensor components. The OME model has the same structure 
except that it contains a flavor (isospin) dependent factor Tj ■ tj. We thus consider the 
following general form of Vqq 

^gg = Y[^i- (^jyA^ij) + CTi • CTj-Ti ■ TjV^r{rij) + SijVririj) + SijTi ■ TjVTr{rij)] , (6) 
with 
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Here, <Tj and Tj are respectively the spin and isospin operators for the ith quark. The radial 
parts of the potentials in Eq. (^) are defined by the momentum space potentials 



^^(0 = / ^Joiqrmq) (8) 



for i = a, ar, and 



V^{r) = J ^32{qr)VM) (9) 

for i = T, Tt. The differences between the OGE model and OME model are in the choices 
of Vi{q) for evaluating Eqs. (H)-®. 



1. OGE Model 
The OGE model is obtained by taking 

K(g) = -^(A.-A,)^F,(g), (10) 
~ Attc\ 1 

^T(g) = -^(A.-A,)^F,(g), (11) 

V^q) = , (12) 
VtM) = , (13) 

where Aj is the color SU(3) generator with 

(Ai ■ \2) = -I , (14) 

for color singlet baryons considered here, and is the quark-gluon coupling constant. In 
Eqs. ([TO|)-(p!l|), we have introduced a form factor of the form 

FM) = (15) 

to regularize the interactions at short distances. This is consistent with the notion that the 
constituent quarks are not point particles within an effective theory. This regularization of 
the gg-potential is essential in obtaining convergent solutions for the bound state problem 
defined by the Hamiltonian hs (Eq. @)). If the potentials are not regularized by form 
factors, the ground state energy is not bound from below. 
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2. OME Model 



Since we only consider baryons with strangeness = 0, the OME model is defined only 
by the exchange of tt and i] mesons. Because r/ is isospin singlet, it only contributes to the 
isospin independent parts of the potential. On the other hand, the exchange of the isovector 
pion leads to isospin dependent terms. The resulting OME model is defined by 

v;.(9) = - f^V i-r^^'>(9) . (18) 



where fuqq is the meson-quark coupling constant, mu denotes the meson mass. Here, as in 
the OGE model, we introduce a form factor 

^M(g) = (20) 

to regularize the potentials at short distances. 

If we assume that the nNN vertex function can be calculated from the nqq interaction 
(as defined later in section III), the vr-quark coupling constant /^gg can be related to the 
Tj'NN coupling constant. Within the naive (Os)'^ configuration for the nucleon, one finds 
that 



f-rrqq _ 3 g^NN 

m7 ~ 5 2M, 



N 



(21) 



where Mjy is the observed mass of the nucleon and we use the empirical value gl^^/An = 14. 
The coupling constant fr,qq can be calculated from /^gg by using the fiavor SU(3) symmetry. 



(22) 



Since there may be SU(3) symmetry breaking mechanisms, the coupling constant /^gg could 
be different from its SU(3) value. We however do not consider this possibility in defining 
the OME model and simply use Eq. (^2]) . 



B. Solution of three-quark bound state problem 

With the Hamiltonian hs (Eq. (|^)) defined above, our first task is to solve the following 
three-body bound state problem in the P = rest frame of the system 
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(23) 



where is the baryon wave function with the label B denoting collectively the spin- 

parity and isospin T ; ms is the mass eigenvalue. We use the diagonalization method 
developed in Ref. ||3^ to solve Eq. ([23|). The basis states for the diagonalization are formed 
from harmonic oscillator wave functions. This choice has two advantages: (1) By using 
appropriate Jacobi coordinates, the center of mass motion can be separated exactly from 
the intrinsic wave functions, (2) The resulting wave functions have the desired S3 symmetry. 
In the diagonalization method, the baryon wave function is expanded as 



(24) 



where the basis wave functions are of the following antisymmetrized form: 



(J) 



■ I'/r 



(25) 



Here a = {a, 



• flavor} stands for the S3 symmetry of each part of the wave function. 



space; '-^^spiin 

l0AfF™) \Xs"^"^) spatial and spin wave functions with the orbital angular momen- 

tum L, quanta of harmonic oscillator A^, and spin S. They couple to give the total angular 
momentum J of the baryon. \ri^^''"'°') and \ip^) are the iso-spin and color wave functions, 
respectively. The color wave function \ip^) is totally antisymmetric. By taking appropriate 
coefficients c^, the basis state | J'^T; i) defined by Eq. (|25D is totally antisymmetric. 

The coefficients a^""^ in Eq. (^) and the mass eigenvalues irtjirrp are obtained from 
diagonalizing the matrix 

H,^, = {rT-i\hB\rT-j). (26) 
In practice diagonalization is performed within a limited number of basis states. Then the 



solution of Eq. (^3]) is a function of the oscillator range parameter h. We treat it as a 
variational parameter and find h by imposing the condition: 



dm 



B 



dh 



. 



(27) 



The basis state is chosen so that the mass eigenvalue rrij^rrp does not change by further 
extension of the basis states. In practice we include the basis states up to llhu. 

For later discussions, we write down here the lowest basis wave functions for A^, A, and 
A^*(S'ii) (simply called A^* from now on) 



I ^^1*5 = 1/2) 
1^^2*5 = 3/2) 



I 2 2' / 
I 2 2' / 



2 2' 



I 1-1. 
I 2 2' 



1) 



2) 



V2 

10^0) 1x^/2) 1^2) IV'^) , 



(1/2) 



,MA\ 



-MA 



MS\ 



MA\ 



(1/2) 



fi') ® IX3^/2)J,/,. \V^}1) + [I^D ® IX^/2)J, \ 1^") , (31) 



,MS\ 



(1/2) 



(1/2) 



(28) 
(29) 

(30) 
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where 0|r^ denotes a state which is totally symmetric with respect to the interchange of any 
pair of quarks, and 0^^^^^^ denote states with mixed symmetries. Similar upper indices are 
also used to specify the symmetry properties of the spin wave function \x) and isospin wave 
function |?7). The spatial wave function can be explicitly written as 

(r, p\<Plo) = [Roo{r)Yo{f) ® i?oo(p)l^o(/))](o) , (32) 

(r, p|0fiS) = _[i?„,(r)Fo(f) ® R,,{p)Y^ip)](^,) , (33) 

{r,p\<pT) = [Ro^{r)Y^{r) ® Roo{p)Yo{p)h) , (34) 

where 

r = n - ra , (35) 

P = r3 — , (36) 

and Rni, Yi are the radial wave function and spherical harmonics respectively. The spin 



and isospin wave functions in Eqs. (p8D-(pl|) can be constructed by using the well known 
procedure. 



III. MESON-BARYON-BARYON VERTEX FUNCTIONS 



Within the constituent quark model, the decay of a baryon (B) into a meson-baryon 
{M'B') state is determined by the matrix element 



(vl/^';M'|i7M(fc)|^^) , 



(37) 



where is a bound state wave function generated from the above structure calculation, and 
HM{k) is an appropriate operator describing how a meson M with a momentum k is emitted 
or absorbed by constituent quarks. In most of the previous works p!6|- [19| , pll , one assumes 
that HM{k) is a one-body operator with the parameters determined phenomenologically by 
fitting some of the partial decay widths listed by PDG. Calculations have also been done 
| 20| by using the ^Pq model for HM{k). 

In this work we assume that HM{k) is a one-body operator which can be derived directly 
from the effective Lagrangian Eq. by taking the nonrelativistic limit of the Feynman 
amplitude ?/p'757'^A;^Mp for the q <-> Mq transition. To be consistent with the nonrelativistic 
treatment of constituent quarks, we keep only the terms up to the order of p/rriq. In 
coordinate space, the resulting q + tc ^ q transition operator is 



H. 



irqq 



k 



2'mn 



[Pi + Pi 



F{k) 



(38) 



where a denotes the z-component of pion isospin and Pi (p'j) is the derivative operator acting 
on the initial (final) baryon wave function; k and a;^ = \Jml + k^ are the momentum and 
energy of pion, respectively. The operator H^igq for the isoscalar rj meson can be obtained 
from Eq. (^) by replacing the label n hj r] and dropping the isospin operator r". 



8 



We note that the operator structure of Eq. 
In Refs. 



is the same as that used in Refs. IIT8|,ET 



18.21 



however the coefficients in front of each term are treated as free parame- 
ters. Here the relative importance between these two terms are fixed by the non-relativistic 
reduction of the effective Lagrangian Eq. (0). To take into account SU(3) breaking, we will 
allow the parameter fr,qq to deviate from its SU(3) value given by Eq. (^2]). Furthermore, we 
also introduce an additional form factor F{k) to account for the effect due to the finite size 
of constituent quarks and mesons. This is consistent with the procedure used above in defin- 
ing the gg-potential Vgq. However, the constituent quark form factor for the interaction Eq. 

could be different from that for the effective gg-potential, since the mesons associated 
with Hj^qq are time-like whereas those associated with Vqq are space-like. The constituent 
quark form factors used in the meson-quark interactions will be introduced in the section 
discussing our results. 



With the operator Eq. 



we find that the M + Bi Bj vertex function can be written 



as 



JM 



(39) 



where Pi, Ji and Tj {pf, Jj and Tf) are the momentum, spin and isospin of the initial (final) 
baryon, respectively. Their 2;-components are denoted by Mj and Tj^ (Mj and Tf^). T, Tz 
and k are the isospin and momentum of meson M (vr or 77); q is the relative momentum of 
the initial meson-baryon system. In the center of mass frame (the rest frame of the final 
baryon Bf), we obviously have the simplification that Pi = —k and q = k. r;g^^. (A;) contains 
the /c-dependence of the vertex function 



ri^B.(^) = 



v/(2vr)32cu^ m 



2J+1 



X -(-.) 



(L100|J0) ( 1 



UJm 



\ Gniq 



m, 



-{Bf\\\3j{lkp)[[Y.j®V\L)®<y\j)r^m) 



(40) 



If the simple wave functions Eqs. (p8|)-(|3lD are used to evaluate Eq.(40), we obtain the 
following analytic expressions for the i:N —>■ N* vertex functions 

(41) 



T^'^ik) = 2T^ik) 



with 



r.(A;) 




v^6V (27r)32cj^ V 



6/ m 



1. 



exp{-y)F{k), 



(42) 



(43) 



where y = 6^/c^/12. Likewise the riN — > A^* vertex functions are found to be 

r??.^(fc) = T%^^{k) = -T,{k). (44) 
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can be obtained from defined in Eq.(43) by replacing tlie label vr by 77. Note that the 
relative importance of the decay vertex functions of and N2 is completely determined 
by the differences in their wave functions given in Eqs. (30)-(31). 



IV. DYNAMICAL MODEL FOR ttN REACTIONS 

With the vertex functions defined by Eq. (^), we follow the procedures of the SL model 
to develop a dynamical model for vrA^ reactions. The starting Hamiltonian is assumed to be 

H = Ho + Hi, (45) 

where the free Hamiltonian takes the following second-quantization form 

Ho = J2 [ dpEBip)biip)bBip) +J2 I dkujM{k)a\^{k)aM{k) ■ (46) 

Here, 6^(65) and a\.j{anj) are the creation (annihilation) operators for the baryons and 
mesons respectively, and 



£b{p) = Jml + p\ (47) 



UMik) = ^Jml + fc2 . (48) 

The baryon mass is generated from the structure calculation described in section II, 
while we use the experimental value for the meson mass rriM- 

The interaction term in Eq. (^SD is written in terms of the vertex functions defined in 
section III 

Hi= I dpdp'dk [(^^'|/7m,,|^''; M)b^s,{p')bB{p)aMik) + h.c] . (49) 

BB'M •' 

The above interaction Hamiltonian is similar to that of the SL model, except that the 



anti-baryon states are absent here. As discussed in Ref. [^, it is a non-trivial many-body 



problem to calculate nN reactions with the use of Hj. To obtain a manageable reaction 
theory, we follow Refs. p9| , p5| and apply the unitary transformation up to the second order 
in Hj to derive an effective Hamiltonian. The essence of the unitary transformation method 
applied in Ref. p9|] is to absorb the unphysical transition B — > M'B' with rriB < rriB' +TnM' 
into non-resonant potentials. The resulting effective Hamiltonian then takes the following 
form 

H,f, = Ho + T + T^ + v , (50) 

where Hq is defined in Eq. (^). The vertex contains only the physical decay process 
B M'B' with rriB > rriM' + iTiB' 

rt = ^ / dkdpdp' (^^'; M'\HM,,\'^'')b^B'{P>iAk')bB{p)0{mB - {niB' + niM')) • (51) 

MBB' 
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where 6{x) = 1(0) for x > 0{x < 0). The non-resonant MB — > M'B' two-body interactions 
are defined by 

v= f dkdk'dpdp\^^'-M'\v\'$^-,M)a\^,{k')aM{k)b^B'ip')bBip)- (52) 

MM' BE' '' 

By translation invariance, the potential matrix element has the following form 

(^^'; M'|{;|^^; M) = 6{p' + k' - p - k){^^'; M>|*^; M). (53) 

To construct the non- resonant interaction v, we are again guided by the SL model. 
We first notice that the low energy Su scattering amplitude can be reproduced very well 
by including the cross nucleon term and p-exchange term. We further notice that the p- 
exchange term in s-wave scattering is equivalent to Weinberg's contact term. Thus for Sn 
scattering considered in this work, it is sufficient to consider the non-resonant mechanisms 
illustrated in Fig. However, we need to extend them to include the transitions to rjN and 
71 A states. 

We write 



V = vt + Vu- (54) 

We derive the u-channel interaction Vu (the first term in the right-hand side of Fig. |I]) by 
using the unitary transformation method presented in Ref. p9|. All transitions between vrA^, 
rjN and vrA states are considered. The resulting matrix elements of w„ are of the following 
form in the center of mass frame: 



Bn 



(55) 



where D{k', k) is given as 



D{k',k) 



1 



+ 



tB{k) + uoM{k) - (e^Jfe + k') + ujM^k) + ujM'{k')) 

1 



(56) 



eB'{k') + ujM'{k') - {tB,Xk + k') + ujM{k) + ujM'{k'))_ 
For ttA^ ttA transition with a nucleon intermediate state, D{k', k) takes a different form 

1 



D{k\k) 



e^ik) + uj^{k) - {e^ik + k') + uo^{k) + uj^{k')) ' 



(57) 



Here B and M denote the baryon and meson states, respectively. An intermediate baryon 
state is denoted by Bn- The allowed intermediate states for each process are listed in Table 
|. The vertex functions in Eq. (^) can be evaluated by using Eqs. (|39|)-(^0D. To obtain 
the partial- wave matrix element from Eq. (|55|), we need to perform the standard angular 
momentum and iso-spin projections. 

For the p-exchange term fj, we assume that it can be replaced by the contact term 
illustrated in Fig. 1 and can be calculated from a contact vrvrgg interaction and the nucleon 
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wave functions generated from the structure calculations. The assumed meson-quark contact 
interaction is of the form 



contact 



(58) 



Taking the matrix element of -/^contact; we obtain 



(vl/^;7r,,|t;i|vl/-;7r,) 



N . 



lewkTkMk) + u^{k'))FN{k - k')Ft{k)Ft{k'). (59) 



Here F^lk) is the iso- vector form factor of the nucleon completely determined by the nucleon 
wave function generated from the structure calculation described in section II. Xt and Ft are 
a phenomenological strength parameter and a constituent quark form factor, respectively. 
These two quantities can be determined by fitting the Su scattering data in the low energy 
region where the N* excitation effects are negligible. They will be given in section VI. 

By using the standard projection operator method |2^, it is straightforward to derive 
from the effective Hamiltonian Eq. ( pUD a calculational framework for vrA^ reactions. The 
transition operator can be written as 



N*,f3 ■ 



(60) 



Here a, (3 denote the meson-baryon states TTN,riN and vrA. A^* are mass eigenstates of 
Eq. (pSl). The first term in Eq. (|60|) is the non-resonant amplitude involving only the non- 
resonant interaction v 



7/3 ' 



(61) 



with 



E - Sb^ (p) - ujal, (k) + ie . 



(62) 



The second term in Eq. (|60D is the resonant term determined by the dressed A^* propagator 
and the dressed vertex functions: 



{E - mN*)6ij - T.ij{E) , 



(63) 
(64) 

(65) 



In Eq. (|63|) , the A^* self-energy is defined by 

s,,(i?) = Er^.,,G;(E)f;^. 



(66) 



The scattering equations defined in Eqs. (|5DD-(|5^) are illustrated in Fig. ^. They are 
solved in the partial-wave representation using the well-known numerical method in mo- 
mentum space. For the Su channel, we consider three meson-baryon channels vrA^, 7]N, vrA 
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and two A^* states. In the vrA channel, we account for the width of the A by modifying the 
propagator Eq. (^) as 



(67) 



where the A self-energy, Sa(co'), is evaluated using the A — ^ vrA^ vertex function determined 



Within the constituent quark model, the nature of the effective gg-interactions has been 
investigated mainly by considering the mass spectrum of nucleon resonances. In this work we 
will apply the reaction model developed in previous sections to further pin down the effective 
gg-interactions by using the ttN scattering data. To see the merits of this approach, it is 
instructive to first consider the simplest case in which the A^, A and N* are described by 
the lowest configurations in the harmonic oscillator basis. The spatial wave functions for 
and A are restricted to s-wave. The A^* states in Su channel are due to Ihu excitation and 
hence there are only two degenerated states |\E'^i'5' = 1/2) and |\E'^2 5' = 3/2). By using 
these simple wave functions given explicitly in Eqs. we are able to obtain analytic 

expressions, which facilitates the understanding of the role of each term in the effective gg- 
interactions, Eq. (^. In particular, the flavor (isospin) structure of the tensor term will be 
shown to be crucial in determining the vrA^ scattering amplitudes. 

With the simple s-wave wave functions Eqs. the mass difference between A^ 

and A is clearly determined solely by the spin-spin interactions of Eq. (^). It is easy to see 



Here (Vi)^ is the matrix element between two gg-states with relative angular momentum 
L = and 1 



The standard notations s and p are used for L = and 1, respectively. For the OGE model 
defined by Eqs. we have Kr = and (K)s > 0. For the OME model defined by 

Eqs. (|16l)-(|l^), we have K- = and (K-r)s < 0. From the signs of the coefficients in Eq. (|68D, 
we can see that both the OME and OGE models can give a positive 6 and can be tuned to 
account for the A-A^ mass splitting. 

For Su states, we diagonalize a 2 x 2 matrix which is obtained by using the wave functions 
l^^r^ = 1/2) and \^^^ S = 3/2) given in Eqs. (|0|)-(|1]) to evaluate Eq. (|2|). Since the 
spin-spin interactions in both the OGE and OME models are of short-range (~ 5-function 
in r-space), we can neglect their matrix elements between p-wave relative wave functions in 
(j)^ and (f)^^. We then find that the difference between the two resulting mass eigenvalues 
has the following analytic form 



in Ref. 



V. RESULTS FROM A SIMPLE MODEL 



6 = niA- rriN = 6(K)s - l2{Var)s- 



(68) 




(69) 



6* 




(70) 



H 



2 
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where A'"^ and denote respectively the lower and higher mass states. The parameter S in 
Eq. (|70|) has already been fixed by the A- mass difference in Eq. (^) . The new parameter 
a is determined by the matrix elements of tensor potentials between two j9-wave relative 
wave functions 

a = -Q{VT)p + 18{VTr)p. (71) 

The resulting wave functions for the two A^* states can be written as 

= cos^|^^i*S = 1/2) + sin^|^^2^ = 3/2), (72) 
= _sine|^^i'S= 1/2) + cos ^1^^2*5 = 3/2). (73) 

The mixing angle 6 also depends on 6 and a 

tan^ = =. (74) 

a + + + 4a2 

The above expressions indicate explicitly how the structure of A is related to that of 
A^* within this simple model. For either of the OGE or OME models, one can adjust their 
coupling parameters to fit the same A-A^ mass difference S. But the difference between 
their tensor potentials will lead to very different a (Eq.(^)), which determines the A^* mass 
splitting 6* (Eq. (^)) and wave functions (Eqs. (0)-(|7|)). We note that the signs of (Vr)p 
evaluated using Eq. (PT| ) for the OGE model and {Vtt)p evaluated using Eq. (P!7| ) for the 
OME model are both positive. It is then clear from Eq. ( fflD that a's for the OGE model 
{Vtt = 0) and for the OME model (Vr = 0) are opposite in sign. Consequently, the phases, 
defined by Eq. ([74D, of the A'£ and A"^ wave functions will be opposite in sign. To see this 
more clearly, we show in Figs. ^§the dependences of 6*/ 6 and mixing coefficient sin^ on the 
parameter a/ 6, which measures the strength of the tensor potential. In the region a/6 < —1 
where the tensor potential resembles that of the OGE model, the A^* mass splitting can be 
smaller than the A^-A mass splitting and the mixing coefficient sin^ is positive. In the 
a/6 > region where gg-potential is close to the OME model, the A^* mass splitting 6* is 
most likely larger than the A^-A mass splitting 6 and the mixing coefficient sin^ becomes 
negative. The results shown in Figs. ^ and ^ clearly indicate that the A^* wave functions 
depend strongly on the fiavor structure of the tensor potential. We can make the OGE, 
OME, or some mixture of OGE and OME models reproduce the same mass splittings 6 and 
6*, but they will yield very different A^* internal wave functions; the difference is particularly 
visible in the relative phases between the l^f^iS* = 1/2) and |\l/^2 5' = 3/2) components. 

We now turn to demonstrating how these different model wave functions can be distin- 
guished by investigating the vrA^ scattering. For this discussion, we neglect the non-resonant 
interaction v of Eq. (|50D and the vrA channel. The irN scattering amplitude in 5*11 channel 
is then determined only by the predicted A^* masses and A^* nN, r]N vertex functions. 
By using the vertex functions defined by Eq. (^) and the A^* wave functions Eqs. (|30D -(pID, 
we obtain in the center of mass frame (p/ = 0,pi = —k, k = q) 

{^^'i^\H^gg\^^;7T) = {-2 COS e + sin e)T^{k) (75) 
(^^^|i7^,,|^^;7r) = (2sin^ + cos^)r^(A;) (76) 
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for the nN N* transition, and 

(^^£|i/,,,|^^;r/) = (cos^ + sin^)r,(A;) (77) 
(*^^|i7^gg|^^;r/) = (-sin^ + cos^)r^(A;) (78) 

for the rjN — ^ A^* transition. Here, r^.^(fc) are defined by Eq. (^31) . The above equations 
clearly show that the couplings of N* states to ttN and rjN continuum are completely 
dictated by 9, which is related to the strength of tensor potential a via Eq. ([7^). This is 
illustrated in Figs. | and |[ We see that in the a/6 < region where the tensor potential is 
close to the flavor independent OGE-type potential (Vy), the lower mass Nl (solid curves) 
decays mainly into the riN channel while the higher mass (dotted curves) favors the decay 
into the ttN channel. The situation is reversed in the a > region where the gg-interaction 
is close to the flavor dependent OME model (Vtt)- 

The results shown in Figs. ^ and |^ suggest that the OGE and OME models will give 
very different vrA^ scattering amplitudes, even their parameters can be adjusted to give the 
same mass splittings 6 = — misi and 6* = rrii^f^ — m7V£- This will be seen in our full 
calculations presented in the next section. 



VI. RESULTS AND DISCUSSIONS 

In this section, we will apply the formulation developed in section IV to explore to what 
extent the commonly used OGE and OME constituent quark models can be consistent with 
the 7iN scattering amplitudes up to 2 GeV. This is obviously a very difficult task since we 
need to consider about 20 partial waves which are known to contain resonance excitations. 
As a start, we will focus on the partial wave. This partial wave involves strong coupling 
between vrA^ and rjN channels and contains two four-star resonances A^(1535) and A^(1650). 
A further complication of this partial wave is that the position of the first resonance is very 
close to the rjN production threshold (Wth = 1485.7 MeV). A detailed study of this channel 
is therefore a very useful first step to get some insights into our approach. 

While we will only focus on the Sn channel, the model must be also consistent with the 
data associated with the well-studied A(1236) resonance. Here we will use the information 



from the SL model ||29| which is consistent with the present formulation and which can 
describe the vrA^ data up to the A excitation. In our interpretation, the bare mass m/^ = 1300 
MeV and bare A — »• ttN form factor determined within the SL model must be reproduced by 
our structure calculations. This is a rather strong constraint on the parameters of the spin- 
spin parts of the residual gg-interactions and the ranges associated with the form factors. 

Another important ingredient in our investigation is the non-resonant interaction v of 
Eq. (|5^). We demand that the constructed non- resonant interactions be consistent with the 
Su amplitude at low energies where the A^* excitation effects are small. This also provided 
a significant constraint in our investigation. The A excitation and low energy vrA^ data were 
not considered in the constituent quark model calculation of vrA^ scattering in Ref. ||30|| . 

In contrast to usual constituent quark model calculations, the determination of the pa- 
rameters in our approach is a highly nonlinear and nonperturbative procedure. For each 
constituent quark model considered, we first carry out extensive structure calculations to 
determine the ranges of its parameters in which the A-A^ mass difference of the SL model 
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can be reproduced. For each possible set of parameters within thus determined ranges, we 
use the predicted wave functions for A^, A, and A^* to calculate various vertex functions 
that arise from the nqq and rjqq interactions. Using these vertex functions, we then calcu- 
late non- resonant potentials. The scattering equations Eqs. (|60|)-(|66|) can then be solved. 
The comparison of the predicted nN amplitudes with the data up to about 2 GeV then 
tells us whether this set of parameters is acceptable. This kind of lengthy structure-reaction 
calculations have to be done many times for each considered constituent quark model until 
the best fit to the data has been obtained. 

In the next few subsections we present our results. 



A. Structure calculations 

The baryon mass eigenvalues and wave functions are obtained by diagonalizing the Hamil- 
tonian Hb defined by Eq. (|^) in the space spanned by the wave functions defined by Eq. 
(^). The variational condition Eq. (pTl) is imposed to determine the oscillator parameter 
h. For all of the models considered in this work, we find that it is necessary to include 
configurations up to IVhu. 

In the one-gluon-exchange (OGE) model defined by Eqs. ([TIl|)-(]r3D, the parameters are 
the vertex cutoff Ag, quark-gluon coupling constant a^,, and the strength of a confinement 
potential. We find that the A-A^ mass splitting depends strongly on the parameters Ag and 
ttg. This is illustrated in Fig. |^ where the nucleon mass is normalized to 940 MeV. We see 
that as as is increased from 0.8 (solid curve) to 1.6 (dot-dashed curve), mA = 1300 MeV 
can be reproduced only when the cutoff Ag is reduced from about 1500 MeV to about 600 
MeV. Similar results are obtained also for other values of the confinement parameter Uc- 

The A^* masses are found to be sensitive to Oc- In Fig. ^ we see that for a wide range of 
tts, the lower mass m7V£ can change by about 200 MeV as Oc is increased from 3 fm~^ (solid 
curve) to 5 fm~^ (dot-dashed curve). However the mass splitting niN^ — rriNi ~ '^^^ MeV is 
less sensitive to Oc, as illustrated in Fig. |^. Since the non-resonant interactions are weak in 
the A^* region, we expect that the parameters must be chosen to yield m7V£ ~ 1500 - 1600 
MeV. From Figs. we see that such values of m^r^ and mA = 1300 MeV can be obtained 
by choosing Ag ~ 1100 MeV, etc ~ 4 fm~^, and ~ 1. 

In the one-meson-exchange (OME) model defined by Eqs. (|I6|)-(|I9]), the parameters are 
vertex cutoffs A^^, A^, and the strength Oc of the confinement potential. Here the meson- 
quark coupling constants f^^gg and /^gg are fixed by Eqs. (PTD-(P2[) using the standard vrA^A^ 
coupling constant g-j^NN and SU(3) symmetry. 

Since the contribution of r^-exchange potential is rather weak, the structure calculations 
are rather insensitive to the range of its cutoff A^. We set A^ = 1000 MeV for simplicity. 
The A-A^ mass splitting is then determined by the parameters A^r and ac alone. In Fig. 



rO| , we see that is rather sensitive to the cutoff A^r. We also find that the Nl mass is 



sensitive to ac, as seen in Fig. |TT]. The dependence of the mass splitting m^v^ — rriN^ on a, 



c 



is shown in Fig. |T2[ We notice here that ttin^ — rriNi ^ 300MeV, which is about a factor 
2 larger than the OGE model value [Fig. ^. This is mainly due to the fact that the OME 
model has a much stronger tensor matrix element, as illustrated in Eq. (^) using the simple 
model. To obtain mA = 1300 MeV and mjv£ ~ 1600 MeV, we need to use A^ ~ 1100 MeV 
and ~ 2 — 3 fm~^. 
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The structure calculations discussed above only identify the possible ranges of the pa- 
rameters for the OGE and OME models. To further pin down the parameters, we now turn 
to discussing our reaction calculations. 



B. Reaction calculations 



As explained in section IV, the first step to perform reaction calculations is to use the 
wave functions from structure calculations to calculate various vertex functions using Eqs. 
(^)- (|40|) . We first notice that all of the predicted A — > vrA^ vertex functions are too hard 
compared with the bare vertex function of the SL model. See the dashed curves in Fig. 



13| . Furthermore, no sensible nN scattering results can be obtained with such hard vertex 
functions. We therefore follow previous works to introduce a constituent quark form factor 
to further regularize the predicted vertex functions. The solid curve in Fig. [13|, which is close 
to the SL model, is obtained by multiplying the dashed curve by the following constituent 
quark form factor 

Another phenomenological aspect in our calculations is to allow the riqq coupling constant 
to vary around its SU(3) value in calculating the rj-BB vertex functions. Namely, we set 



(80) 



where Xj^qq is a phenomenological parameter that is allowed to vary along with the param- 
eters ko and [Eq. (|7|)] in fitting the nN amplitudes. 

The next step is to fix the non-resonant potential Vt defined by Eq. (^) in the low energy 
region where the N* excitation effects are negligible. We find that the data up to ~ 1300 
MeV can be described well if we take Xt = 0.7 and set Ft{k) as a dipole form with a 500 MeV 
cutoff mass. A monopole constituent quark form factor with 1 GeV cutoff is also included 
to soften the vertex functions of the u-channel interaction (the first term of Fig. |1|). 

We have found that both the OGE and OME models, as defined in this work, can give a 
good description of the Su amplitudes up to ly < about 1500 MeV. They however can not 
describe the data at higher energies. The best results we have obtained are shown in Figs. 
I4|-p!^. We see that the OGE model can reproduce the rapid change in phase at ~ 1500 



MeV, while the OME model fails completely. The resulting parameters are listed in Tables 
|I-III for the structure calculations and in Table IV for the calculations of A^* ttN, rjN, 



IT A form factors. 

To understand the differences between the OGE and OME models, we show in Figs. 
|T6|-|1^ the vertex functions calculated using the best-fit parameters listed in Tables ^ and 
p| . Generally the spin-fiavor structure of gg-interaction strongly infiuences the predicted 
vertex functions. For the A^* ttN vertex functions (Fig. ^), and A"^ have almost 
the same strength in the OGE model, while decays more strongly than A^^ in the OME 
model. This can be understood by simply considering the two main Ihu components of 
A^* given in the first two rows of Table 0. The main difference between the wave functions 
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of the OGE and OME models is in the relative sign between the S = 1/2 and S* = 3/2 
components. This is due to the flavor structure of tensor potential as illustrated in Fig. ^ in 
section V. This difference in the structure of A^* wave functions also plays an important role 
for the A^* — >■ rjN strength, as seen in Fig. 0. In the OGE model the strength of A'^^ — rjN 
is stronger than —>■ rjN, while the situation is opposite in the OME model. For the 
A^* —>■ ttA vertex functions compared in Fig. |18|, we again see a large difference between the 
two models. The OME model predicts a very weak strength for the higher mass N^. 

In a high energy region W > 1500 MeV, we find that the predicted vrA^ scattering 
amplitudes are dominated by the resonant term of Eq. (0). Thus the energy dependence 
of scattering amplitudes in this region can be clearly understood by examining the A^* 
propagator D{E) defined by Eq. (|^). Here we see that the coupling of the A^* states to 
nN, r]N and ttA continuum can shift their masses by S(-E'). The poles of D{E) are the 
resonance positions. These can be obtained by diagonalizing D{E) in the space A'^^ © A'"^. 
The real parts of the resulting poles include the mass shifts due to the couplings between the 
two A^* states via meson-baryon continuum, and the imaginary parts are the widths of the 
resonances. It is important to note here that these effects due to meson-baryon continuum 
are not included in the structure calculations described in section II or any of the existing 
constituent quark model calculations. In Fig. O we show the resulting mass eigenvalues 
of D{E). The intersections between the mass eigenvalues and the dotted lines representing 
Re(i?*) = W are the resonance positions. In our fit, X^gg and the constituent quark form 
factors are adjusted for each model so that the resonance position for A'^^ lies around 1535 
MeV, the PDG value. However, the predicted position for depend very much on the 
model. In the OME model the resonance position of is ~ 1900 MeV which is clearly too 
high compared with the PDG value, 1650 MeV. Even in the OGE model it is still too high 
~ 1800 MeV. This is the main difficulty we have in fitting the data at high energies. If we 
choose parameters that fit a lower resonance position for N^, becomes too light and the 
data below W < 1500 MeV can not be fitted at all. 

In addition to the resonance positions, the OGE and OME models differ also in the 
energy dependence of the imaginary parts of the self-energy S(-E). This is illustrated in Fig. 



20| . The strong energy dependence of the imaginary part for A'^^ of the OGE model is due 
to its stronger coupling with the rjN channel, as seen in Fig. 0(a). This leads to a very 
strong energy dependence due to the opening of rjN threshold; see the solid curve of Fig. 
pO|(a). The strong coupling of to the riN channel is essential in reproducing the rapid 
change in phase shift around W = 1500 MeV. This feature can not be generated by the 
OME model within which the rjN channel is mainly coupled to the higher mass A^^, as can 
be seen in Fig. P^(b) and |2D|(b). This is why the OGE model is better than the OME model 
in reproducing the data up to about 1550 MeV. 



C. Phenomenological model 

The above results are qualitatively consistent with the results of the simple model given 
in section V. Our findings are perhaps consistent with a recent phenomenological study 
T7| of negative parity nucleon resonances. It was also found there that the effective qq- 



interaction is not simply given by the OME mechanism. As an attempt to improve the fit 
to the 5*11 amplitude, we have also explored the mixture of OGE and OME models. It turns 
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out that such a hybrid model also fails, mainly due to the very disruptive tensor component 
of the OME model in determining the phases of wave functions, as discussed in section V. 

We therefore turn to investigating a purely phenomenological model. We first observe 
that the reason why the OGE model is better than the OME model in reproducing the main 
features of the energy- dependence of the data is that its tensor force yields correct relative 
phases between the two A^* wave functions. Consequently, a successful phenomenological 
model should have a flavor independent tensor force similar to that of the OGE model. 
Second, we observe that the difficulty of the OGE model in reproducing the data at higher 
energies is mainly due to its very large mass splitting mjy^ — rriN*. As discussed in section 
V within the simple model, this problem can not be fixed within the OGE model unless the 
A^-A mass splitting is not constrained by the SL model. This difficulty can be avoided if 
we have a short range flavor- dependent spin-spin interaction like that of the OME model. 
Qualitatively speaking, the data of A excitation and Sn nN scattering seem to favor a tensor 
term due to one-gluon-exchange and a spin-spin interaction due to one-meson-exchange. 

The above considerations have guided us to explore many phenomenological models. For 
example, we have found that the vrA^ 6*11 amplitudes can be much better described by the 
following phenomenological model 

K = 0, (81) 
y<^r = ^^Ph(^), (82) 

Vt = ^FM, (83) 

Vtt = 0. (84) 

The resulting parameters are listed in Table |V1[ It is interesting to flrst compare the A^* 
wave functions of this model with those of the OGE model. We see in Table that the 
relative phases between the S = 1/2 and S = 3/2 components in the A^* wave functions 
are the same in the two models. We therefore expect that the A^* — vrA^, rjN, vr A vertex 
functions also must be qualitatively very similar. However, we see from Table that the 
phenomenological model apparently has a stronger tensor potential in mixing the S = 1/2 
and 5 = 3/2 components. As discussed in section V using the simple model, we therefore 
expect that the differences between the and A"^ vertex functions must be larger than 



that of the OGE model. This is exactly what we see by comparing the results in Fig. ^ 
and Figs. |g(a), |T3(a), and 0(a). 

According to Tables y and |VI|, rriN^ = 1719.5 MeV for the phenomenological model is 
lower than the value 1772.4 MeV of the OGE model. This about 50 MeV shift of mjv|^ is 
also crucial in improving the flt at higher energies. 

The results from the phenomenological model are compared with the nN data in Figs. 



22|-p^. The predicted vrA^ Sn phase shifts and Su amplitudes become closer to the data than 
the OGE model. Moreover the rapid energy- dependent behavior of the phase shift around 
W ~ 1500 MeV is reproduced. This behavior is closely related to the sharp structure in the 
width of the in the solid curve of Fig. ^(b). We also observe that Fig. ^(b) indicates 
that below the 77 A^ threshold ~ 1490 MeV, the imaginary part of D(E) for A^^ is smaller 
than that of N^. The situation is opposite in the OME model (Fig. pO|(a)). This change in 
the widths is also instrumental in getting better results at energies near riN threshold. 
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The improvement due to the A^* mass sphtting is clearly reflected in the nN scattering 
amplitudes in Fig. 23. Compared with the OGE model results (solid curves in Fig. |l^), the 
peak positions are much better reproduced. However, there are still significant discrepancies 
near the second peak. This can be understood from the pole positions of T-matrix, displayed 
in Fig. The position of the second resonance is at about 1700 MeV, still 50MeV higher 
than the PDG value, 1650 MeV. On the other hand, the first resonance position is very close 
to the PDG value, 1535 MeV. 

Finally the total cross sections of the ti" + p ^ rj + n reaction predicted by using our 
three models (OME, OGE and phenomenological) are compared with data in Fig. ^ Near 
the rjN production threshold, where is expected to dominate the cross section, both 
the OGE and phenomenological models explain the data well, while the results from the 
OME model are too small due to the weak coupling of its A'^^ to the r]N channel. The 
phenomenological model seems to give the best description of both the ttN amplitudes and 
the 7] production cross sections. The discrepancy with the data at higher energies in Fig. ^ 
is mainly due to the neglect of non-Sn partial waves. 

Our results in this section suggest that the residual gg-interactions within the constituent 
quark model are much more complicated than the conventional OGE and OME models. 
Within the effective theory defined by the Lagrangian Eq. (|^), higher order exchanges of 
mesons and gluons must be considered. The phenomenological model we have obtained is 
very suggestive. It remains to be seen whether this model is also consistent with the data 
of other partial waves. 



VII. CONCLUSIONS 

We have developed a dynamical approach to predict ttN scattering amplitudes starting 
with the constituent quark models. This is an extension of the dynamical model developed 



in Ref. |]2^ to account for the multi-channel and multi-resonance cases. In this exploratory 
investigation, we focus on the ttN amplitude in the Sn channel and only consider the most 
frequently used nonrelativistic constituent quark models based on either the one-gluon- 
exchange (OGE) or the one-meson-exchange (OME) quark-quark residual interactions. 

The first step of our calculations is to choose appropriate parameters of the considered 
constituent quark models to reproduce the bare parameters associated with the A within the 



ttA^ model developed in Ref. |^9[. Here, we apply a variational method to solve the three- 
quark bound state problem. The resulting wave functions are then used to calculate the 
N* —>■ ttN, 7]N, ttA vertex functions by assuming that the tt and r] mesons couple directly 
to quarks. These vertex functions and the predicted baryon bare masses then define a 
Hamiltonian for irN reactions. We apply the unitary transformation method of Ref. p9[ to 



solve the vrA^ scattering problem. The final parameters of the considered constituent quark 
models are determined by fitting the nN scattering data. 

We have found that both the OGE and OME models can reproduce the Sn scattering 
amplitudes only up to about W = 1500 MeV. The OGE model is better in reproducing the 
rapid change in ttN phase near the t] production threshold, as shown in Fig. ^ However, 
both models fail to describe the data at higher energies. The dynamical origins of the 
difficulties are found to be due to the sensitivities of the predicted A^* —>■ ttN, rjN vertex 
functions to the structure of the assumed residual quark-quark interactions. In particular, it 
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is found that the flavor dependent iji-Tj) tensor component of the OME model does not seem 
to be favored by the data. On the other hand, the OGE model can describe the data better 
if its spin-spin interaction includes a flavor dependent factor. To illustrate this, we have 
shown that the data can be reasonably well described (Figs. ^^ - |23|) by a phenomenological 
model which has a spin-spin interaction from the OME model and a tensor interaction from 
the OGE model. 

In conclusion, our results indicate that the residual quark-quark interactions within the 
nonrelativistic constituent quark model could be much more complicated than the simple 
OGE and OME mechanisms. In the future, we need to consider relativistic effects and the 
residual quark-quark interactions due to multi-gluon and/or multi- meson exchanges. In ad- 
dition, we need to investigate the two-body effects on the calculations of A^* nN, t]N, tt A 
vertex functions. Our investigations in these directions will be published elsewhere. 
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TABLES 



TABLE I. Baryon states included in the u-channel non-resonant potential. 
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TABLE II. Parameters and results of masses in the OGE model. 

Oc [fm-^] Os Ag [MeV] M(Af£) M(iV^) 

lo LO 1087 1593.5 1772.4 



TABLE III. Parameters and results of masses in the OME model. 

«c K \ M{Nl) M{N}j) 

2.5 1139 1000 1565.0 1885.8 



TABLE IV. Parameters of Fermi type form factors. 
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TABLE V. Coefficients of Ifiu; configurations. 
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TABLE VI. Parameters and results of masses in the phenomenological model. 

oc Oar or Aph M (Nl) MjNlj) 

2.5 0.427 1.60 1200 1580.6 1719.5 
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FIGURES 
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FIG. 1. Graphical representation of the non-resonant meson-baryon interaction. 
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FIG. 2. Graphical representation of the scattering equations defined by Eqs. (|60|) -(|66|). 
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FIG. 3. Mass splitting 5*(Eq.(70)) of N* as a function of the matrix element a(Eq.(71) of the 
tensor potential. Both are in units of A- AT mass difference 6. 
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FIG. 4. Mixing coefficient of the N* eigenfunctions (Eqs.(72)-(73)) as a function of the matrix 
element a(Eq.(71)) of the tensor potential. 
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FIG. 5. N* rjN strength as a function of the matrix element a(Eq.(71)) of the tensor 
potential. The solid and dashed curves are the strengths for A^£ and A^^, respectively. 
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FIG. 6. Same as Fig. |5| except for the vrA^ transition. 




900 



400 



800 1200 1600 2000 
Ag [MeV] 



FIG. 7. Mass of A as a function of the cutoff Ag in the OGE model. The solid, dashed and 
dot-dashed curves are results with = 0.8,1.0 and 1.6, respectively. Here the nucleon mass is 
normalized to 940 MeV and = 4 fm~^ is used. 
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FIG. 8. Mass of N'l as a function of in the OGE model. The sohd, dashed and dot-dashed 
curves are results with Oc = 3, 4 and 5 fm~^. For each value of ac and a^, the cutoff parameter Ag 
is determined by requiring the mass of A to be 1300 MeV. 
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FIG. 9. Mass splitting of N*'s in the OGE model. See caption of Fig. 
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FIG. 10. Mass of A as a function of A^r in the OME model. The sohd, dashed and dot-dashed 
curves are the results with = 2, 3 and 4 fm~^ and A„ = 1 GeV. 
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FIG. 11. Mass of iV£ as a function of in the OME model. The cutoff for each Oc is 
determined so that the mass of A is 1300 MeV. 
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FIG. 12. Mass splitting of N* in the OME model. See caption of Fig. 11 
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FIG. 13. A ^ ttAT vertex function in the OGE model (a) and the OME model (b). The 
solid and dashed curves are the vertex functions with and without the quark form factor 
respectively. The dot-dashed curve is the vertex function of the SL model. 
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FIG. 14. Phase shifts of the vrA^ scattering in Sn channel. The sohd and dashed curves are 
the results from the OGE and OME models, respectively. Open circles are the data of VPI SP98 
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FIG. 15. Same as Fig. 14 except for the real parts (a) and the imaginary parts (b) of the vrA'^ 
scattering ampUtudes in Sn channel. 
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FIG. 18. Same as Fig. 16 except for the A^* — > vrA vertex functions. 
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FIG. 19. Real parts of the eigenvalues of D{E){Eq.{63)) calculated from the OGE model (a) 
and the OME model (b). The solid and dashed curves correspond to the masses for and N^, 
respectively. The dotted line represents Re{E*) = Wcm- 
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FIG. 20. Same as Fig. ^ except for the imaginary parts of the eigenvalues of D{E). 
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FIG. 21. AT* ttN (a), A?^* — 77A/" (b), and iV* — ttA (c) vertex functions in the phenomeno- 
logical model. The solid and dashed curves are for and N^, respectively. 
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FIG. 22. Phase shifts of the vrA^ scattering in Su channel calculated from the phenomenological 
model. Open circles are data from VPI SP98 
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FIG. 23. Same as Fig. 22 except for the real parts (a) and imaginary parts (b) of the vrA^ 
scattering amplitudes of in 
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FIG. 24. Real parts (a) and imaginary parts (b) of the eigenvalues of D{E)(Eq.{63)) calculated 
from the phenomenological model. The solid and dashed curves correspond to the masses for 
and N^, respectively. 
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FIG. 25. Total cross section ofvr + p ^ rj + n reaction. The dashed, dot-dashed and sohd 



curves are results of the OME, OGE and phenomenological models. Data are taken from Ref. |36|. 
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